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Abstract 

This is a brief overview of RMT applications to quantum or wave chaotic reso- 
nance scattering with an emphasis on non-perturbative methods and results. 

34.1 Introduction 

Describing bound states corresponding to discrete energy levels in closed quan- 
tum systems, one usually addresses properties of Hermitian random matrices 
(Hamiltonians) H. Experimentally, however, one often encounters the phe- 
nomenon of chaotic scattering of quantum waves (or their classical analogues) 
in open systems |Kuh05a] . The most salient feature of open systems is the set 
of resonances which are quasibound states embedded into the continuum. The 
resonances manifest themselves via fluctuating structures in scattering observ- 
ables. A natural way to address them is via an energy-dependent scattering 
matrix, S(E), which relates the amplitudes of incoming and outgoing waves. 
In such an approach, the resonances correspond to the poles of S(E) located 
(as required by the causality condition) in the lower half-plane of the complex 
energy plane |Nus72]. 
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The correspondence becomes explicit in the framework of the Hamiltonian 
approach [Ver85]. It expresses the resonance part of the S-matrix via the 
Wigner reaction matrix K(E) (which is a MxM Hermitian matrix) as follows 

3(E) = \ ~ , K(E) = lV^E-Hy i V, (34.1.1) 

where V is an NxM matrix of energy-independent coupling amplitudes be- 
tween N levels and M scattering channels. By purely algebraic transforma- 
tions, the above expression can be equivalently represented in terms of the 
NxN effective non-Hermitian Hamiltonian T~L e s of the open system: 

S(E) = l-iVi * V, n cS = H-iVVK (34.1.2) 

hi - HeS 

The anti-Hermitian part of H e g has a factorized structure which is a direct 
consequence of the unitarity condition (flux conservation): S'(E)S(E) = 1 for 
real E. It describes a coupling of the bound states via the common decay 
channels, thus converting the parential levels into N complex resonances £ n = 
E n — |T n characterized by energies E n and widths T n > 0. The eigenvectors 
of 1~L c ff define the corresponding bi-orthogonal resonance states (quasimodes). 
They appear via the residues of the 5-matrix at the pole positions. 

Universal properties of both resonances and resonance states in the chaotic 
regime are analysed by replacing the actual non-Hermitian Hamiltonian with 
an RMT ensemble of appropriate symmetryQ In the limit N — > oo spectral 
fluctuations on the scale of the mean level spacing A turn out to be universal, 
i.e. independent of microscopical details such as the particular form of the 
distribution of H or the energy dependence of A. Similarly, the results are also 
independent of particular statistical assumptions on coupling amplitudes V£ as 
long as M <C N |Leh95a| . The amplitudes may be therefore equivalently chosen 
fixed [Ver85] or random [Sok89] variables and enter final expressions only via 
the so-called transmission coefficients 

Tc = 1 - \S CC \ 2 = , AK \ , k c = (34.1.3) 
1 1 (1 + k c ) 2 ' 2iVA v ' 

where S stands for the average (or optical) S matrix. The set of T c 's is assumed 
to be the only input parameters in the theory characterising the degree of system 
openness. The main advantage of such an approach is that it treats on equal 
footing both the spectral and scattering characteristics of open chaotic systems 
and is flexible enough to incorporate other imperfections of the system, e.g., 
disorder due to random scatterers and/or irreversible losses |Fyo05| . 



Conventionally, H is taken from the GOE or GUE labeled by the Dyson's index /3 = 1 
or P = 2, respectively, according to time-reversal symmetry being preserved or broken in the 
system (GSE, j3 = 4, is to be taken when spin becomes important). 
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In fact, the description of wave scattering outlined above can be looked 
at as an integral part of the general theory of linear dynamic open systems 
in terms of the input-output approach. These ideas and relations were devel- 
oped in system theory and engineering mathematics many years ago, going 
back to the pioneering works by M. Livsic |Liv73| . A brief description of main 
constructions and interpretations of the linear open systems approach, in par- 
ticular, a short derivation of Eq. (134.1.2p . can be found in |FyoOO| . In what 
follows, we overview briefly selected topics on universal statistics of resonances 
and scattering observables, focusing mainly on theoretical results obtained via 
non-perturbative methods starting from mid-nineties. For a more detailed dis- 
cussion, including applications, experimental verifications as well as references, 
we refer to the special issue of J. Phys. A 38, Number 49 (2005) on 'Trends in 
Quantum Chaotic Scattering' as well as to the recent reviews |Bee97t tGuh98, 
IMel99j lAlhOOl IMitlOj . Finally, it is worth mentioning that many methods and 
results discussed below can be extended beyond the universal RMT regime and 
the use of more ramified random matrix ensembles (e.g. of banded type) al- 
lows one to take into account effects of the Anderson localisation. Discussion 
of those developments goes, however, beyond the remit of the present review, 
and an interested reader may find more information and further references in 
|Kot02l [EVo03b| IMen05l IKot05l lOssOSl |Fyo05t IWei06l IMir061 lMon09j . 



34.2 Statistics at the fixed energy 



Statistical properties of scattering observables considered at the given fixed 
value of the scattering energy E can be inferred from the corresponding distri- 
bution of S = S(E). The later is known to be uniquely parameterized by the av- 
erage S"-matrix and is distributed according to Poisson's kernel |Hua63[lMel85] : 



Ps(S) 



1 



det[l 



S ] S] 



det[l 



■ S ] S] 2 



(PM+2-P)/2 



(34.2.1) 



where Vp is a normalization constant. Although this expression follows from 
Eq. (|34.1.2p in the limit N — > oo, as shown by Brouwer [Bro95| . Pg(S) can be 
derived starting from statistical assumptions imposed on S without referencing 
to the Hamiltonian at all, as it was initially done by Mello et al. [Mel85j . For a 
detailed account of the theory see the recent review [Mel99] and book |Mel04| . 



34.2.1 Maximum entropy approach 

The starting point of that information-theoretical description is an information 
entropy S = — f dfi(S) Pg(S) In Pg(S) associated with the probability distribu- 
tion P-g(S). The integration here is over the invariant (Haar) measure dfi(S) 
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which satisfies the symmetry constraints imposed on S: unitary S should be 
chosen symmetric for (3 = 1 or self-dual for /3 = 4. Due to the causality condi- 
tion, S(E) is analytic in the upper half of the complex energy plane, implying 
S k = S k for the positive moments, where the average is assumed to be per- 
formed over an energy interval. One can further assume that S(E), for real E, 
is a stationary random-matrix function satisfying the condition of ergodicity, i.e. 
ensemble averages are equal to spectral averages, the later being determined by 
the (given) optical scattering matrix S. These two conditions together yield the 
so-called analyticity- ergodicity requirement J dfj,(S)S Pg(S) = S , k = 1, 2, . . .. 
The sought distribution is then assigned to an "equal a priory" distribution 
of S, i.e. the one which maximizes the entropy S, subject to the above con- 
straints. The answer turns out to be exactly given by Eq. (|34,2.ip [Mel85], with 
V/3 = J dfj,(S) being the volume of the matrix space. Then for any analytic func- 
tion f(S) its average is given by the formula f dfi(S) f (S)P-g(S) = /(S). Thus 
the system-specific details are irrelevant, except for the optical S-matrix. 

The important case of S = (ideal, or perfect coupling) plays a special role, 
being physically realised when prompt or direct scattering processes described 
by S are absent. The case corresponds to Dyson's circular ensembles of random 
unitary matrices S distributed uniformly with respect to the invariant measure 
(note Pq(S) = constant). Thus statistical averaging amounts to the integration 
over unitary group, various methods being developed for this purpose [Mel85, 
Bro96 1 . The general situation of nonzero S is much harder for analytical work. 
In this case, however, one can use the following formal construction 

S = (ti)- 1 ^ - S)(l - tfS)-H\ , (34.2.2) 

where So is drawn from the corresponding circular ensemble and the matrices t\ 
and t[ belong to a 2Mx2M unitary matrix Si = (J 1 r )) of the same symmetry. 
Evaluating the Jacobian of the transformation (|34.2.2p , one finds |Hua6 3 , Fri85 
that the distribution of S is exactly given by the Poisson's kernel (|34.2.ip pro- 
vided that we identify r\ = S. Physically, this can be understood by subdivid- 
ing the total scattering process into the stage of prompt (direct) response Si 
occurring prior to the stage of equilibrated (slow) response So |Bro96l IMel99j , 
Eq. (|34,2.2p being a result of the composition of the two processes. In this way 
one can map the problem with direct processes to that without and several 
useful applications were considered in |Gop98 , ISavOlj Fyo05 , Gop08 1 . 



In problems of quantum transport it is often convenient to group scattering 
channels into N\ 'left' and N2 'right' channels which correspond to the propa- 
gating modes in the two leads attached to the conductor (with M = N\ + N2 
channels in total). In this way the matrix S acquires the block structure 
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where r (r') is a AixAi (A^xA^) reflection matrix in the left (right) lead and 
t (t') is a rectangular matrix of transmission from left to right (right to left). 
Unitarity ensures that the four Hermitian matrices tt\ t't'\ l — rr^ and 1 — r'r'* 
have the same set of n = mu^A^A^) nonzero eigenvalues {Tk G [0,1]}. In 
terms of these so-called transmission eigenvalues the S matrix can be rewritten 
by employing the singular value decomposition for t = vtv! and if = ut t v' 
yielding the second equality in (134. 2. 3p . Here t is a JViX A2 matrix with all 
elements zero except Tkk = VTk, k = l,...,ra, u and v! (v and v') being 
N±xNi (N2XN2) unitary matrices (u*vf = v*v' = 1 for /3 = 1,4). For S drawn 
from one of Dyson's circular ensembles those unitary matrices are uniformly 
distributed (with respect to the Haar measure), whereas the joint probability 
density function (JPDF) of the transmission eigenvalues reads [Bar94l [Jal94] 

1 n R 

MiW = TrWfflf U 7 ?- 1 , c = k\Ni~ N2\ + 1) • (34.2.4) 

Here and below A({xj}) = ni<?( 

) denotes the Vandermonde determi- 
nant, A//3 being the normalization constant. Since the Landauer-Biittiker scat- 
tering formalism allows one to express various transport observables in terms of 
Tj's |Bla00| . Eq. (|34.2.4[) serves as the starting point for performing statistical 
analysis of quantum transport in chaotic systems, see |Bee 97] for a review. 

34.2.2 Quantum transport and the Selberg integral 

The conductance g = tr(tt') = Ya=i Ti and the (zero-frequency) shot-noise 
power p = Ya=i 710- ~ Ti) are important examples of a linear statistic on trans- 
mission eigenvalues. Their mean values are provided by averaging over t or 
by making use of the mean eigenvalue density (derived recently in [Viv08]). 
However, the higher moments involve products of 7}'s and require to integrate 
over the full JPDF. This task can be accomplished by recognizing a profound 
connection of (|34.2.4p to the celebrated Selberg integral [Sav06aj which is a 
multidimensional generalization of Euler's beta-function^ The specific struc- 
ture of its integrand yields a closed set of algebraic relations for the moments 
(7i Al ' ' ' Tm m ) °f given order r = ^ Aj with integer parts Aj > |Sav08| . 

As a particular example, we consider the relations for the second order mo- 
ments: [a + 2 + p(n - l)}(7?) = [a + I + (3{n - l)](Ti) - f (n - l){Ti%) and 
(71 • • • T m ) = njLi a+ i+p(2n-j-i)/2 at m = 2. These provide us straightfor- 
wardly with the well-known exact results for both the average and variance of 

2 / 1 dT 1 .../ i dr n |A(T)i 2 '=nr= 1 T; o - 1 (i-T i ) 1 '- 1 = n;^ 1 m being 

the gamma function, gives the definition of Selberg integral valid at integer n > 1, complex 
a and b with positive real parts, and complex c with Rec > — min[A, Re ^rrj, Re — rr]. The 
constant Mp is given by this expression at a = a, b = 1 and c= %. Chapter 17 of Mehta's 
book gives an introduction into the Selberg integral theory, see [For08| for its current status. 
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the conductance (first obtained by different methods in [Bar941 lJ al94[ [Bro96j) 
as well as with the average shot-noise power (hence, Fano factor |Sav06a| : 

/ \ N,N 2 (^-1 + 1)^-1 + 1) 

i9) = mTT3|' {P) = i9) (M-2 + |)(M-l + i) - (342 - 5) 

(We note also an interesting exact relation var(g) = ^(NiN 2 )^ 1 (g) (p) .) The 
approach has been developed further to study full counting statistics of charge 
transfer [Nov07| as well as to obtain explicit expressions for the skewness and 
kurtosis of the charge and conductance distributions, and for the shot-noise 
variance jSav08] . With symmetry index /3 entering the Selberg integral as a 
continuous parameter, the method allows one to treat all the three ensembles 
on equal footing and provides a powerful non-perturbative alternative to orthog- 
onal polynomial [Ara981 IViv08| or group integration approaches [B ro9 6| IBul06| . 

It is further possible to combine Selberg integral with the theory of symmet- 
ric functions and to develop a systematic approach for computing the moments 
of the conductance and shot-noise (including their joint moments) of arbitrary 
order and at any number of open channels [Kho09| . The method is based on ex- 
panding powers of g or p (or other linear statistics) in Schur functions s\. These 
functions are symmetric polynomials s\(T) = det{7^ A:,+ra ^}™j = i/jiet{T™ ~*}fj= 
in n transmission eigenvalues 71, ... ,T n indexed by partitions A a In the group 
representation theory the Schur functions are the irreducible characters of the 
unitary group and hence are orthogonal. This orthogonality can be exploited 
to determine the "Fourier coefficients" in the Schur function expansion by in- 
tegration over the unitary group. The general form of this expansion reads: 

n i +oo x 

Il( E a i T A = I>A(a>A(T), c A (o) = det{a Xk -k+i}l l=r (342 - 6) 

i=l S'=-oo ' X 

The summation here is over all partitions A of length n or less, including empty 
partition for which s\ = 1. The Schur functions can be then averaged over the 
JPDF (|34.2.4p with the help of integration formulas due to |Hua63j . One gets 

^(A j + Nx-j)] (A j+ N 2 -j)\ (M-Z(A)-i)! 



; , (m-jy. (N 2 -j)\ (\ j+M -i(\)-j)\ 



l<i<j<l{X) J 



3 A partition is a finite sequence A = (Ai, A2, . . . , A m ) of non-negative integers (called parts) 
in decreasing order Ai > A2 > . . . > A m > 0. The weight of a partition, |A|, is the sum of its 
parts, |A| = \j, and the length, Z(A), is the number of its non-zero parts. 
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where the coefficient c\ = ^^'fAf'^^ J ~— - has been introduced, and 

, . ^ (A J +iV 1 -j)!(A J +iV 2 - J )! (M-j)! 

(SA> ^ 2 " CA /J (M-i)! (AW)! (A,+Af- J )!' (34 - 2 - 8) 

in the cases of orthogonal (/3 = 1) and unitary {(3 = 2) symmetry, respectively. 

As an application of the method we consider the moments of eg + p. This 
quantity has a physical meaning of the total noise including both thermal and 
shot-noise contributions [BlaOO] . From expansion (jMjZBD one finds [Kho09| the 
r-th moment of the total noise given by 



r 



((eg+py} = r\Y J (-l) m a + ey- m Y J A,m(«A>, (34-2.9) 

m=0 \X\= r +m 

where the second sum is over all partitions of r+m. Here the expansion coeffi- 
cient fx >m = det{ fc .j^,_jj ) _ 3 -_2fc .); }) where the summation indices ki run over 
all integers from to m subject to the constraint k\ + . . .+fcjm = m. In the limit 
e — > oo Eq. (I34.2.9|) yields all the conductance moments (g r ) = r\ J2\x\=r c a(sa) 
(note that /a,o = cx introduced above) first obtained in |Nov08| at j3 = 2 (see 
[Osi08| for an alternative treatment). In the opposite case of e = one gets 
all the moments of shot-noise, hence all cumulants by well-known recursion re- 
lations. Apart from producing lower order cumulants analytically, the method 
can be used for computing higher order cumulants symbolically by employing 
computer algebra packages. In particular, one can perform an asymptotic anal- 
ysis of the exact expressions in the limit N\ 2 3> 1 to make predictions for the 
leading order term in the i expansion of higher-order cumulants [Kho09j . 

As concerns the corresponding distribution functions, explicit expressions 
for the conductance distribution Pg^\g) are available in particular cases of 
n = 1, 2 [Mel99| . At N 2 = K > N x = 1, Eq. (134.2.41) readily gives P^ =1 (g) = 
(pK/2)gP K / 2 - 1 , < g < 1. In the case of N\ = 2 and arbitrary N2 = K > 2, 
one can actually find an exact result at any positive integer (3 |KhoQ9j which 
reads: P^ 2 (g) = Kg^-\X X - (-l)W*-Mx 2 @(g- 1) £?=o {^B X ^{K- 
1) + j, 1 — /3K)] for < g < 2 and zero otherwise. Here B z (a, b) stands for the 

incomplete beta- function, X\ = ^^^/^r^xf^ 2 ^ an< ^ ^ 2 = T(p)r[p(K-i)%\ • 
As both N\ and ^2 grow, the distribution of the conductance or, generally, of 
any linear statistic rapidly approaches a Gaussian distribution [Bee93]. With 
exact expressions for higher-order cumulants in hand, one can obtain next order 
corrections to the Gaussian law by making use of the Edgeworth expansion 
|Bli98] . Such approximations were shown to be fairly accurate in the bulk even 
for small channel numbers. Actually, one can derive an explicit representation 
for the distribution functions in terms of Pfaffians. It works in the whole density 
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support including the spectral edges where the distributions have a power-law 
dependence. Such asymptotics can be investigated exactly for any f3 and Ni t %, 
including powers and corresponding pre- factors, see [Kho09] for details. 



34.3 Correlation properties 

The maximum-entropy approach proved to be a success for extracting many 
quantities important for studying electronic transport in mesoscopic systems, 
see Chap. 35. At the same time, correlation properties of the S-matrix at 
different values of energy E as well as other spectral characteristics of open 
systems related to the resonances turn out to be inaccessible in the framework 
of the maximum-entropy approach because of the single-energy nature of the 
latter. To address such quantities one has to employ the Hamiltonian approach 
based on (j34.1.2|) . Being supplemented with the powerful supersymmetry tech- 
nique of ensemble averaging, this way resulted in advances in calculating and 
studying two-point correlation functions of S'-matrix entries and many other 
characteristics, as discussed below. However, exact analytical treatment of the 
higher-order correlation functions (needed, e.g., in the theory of Ericsson fluc- 
tuations, see Chap. 2) remains yet an important outstanding problem. 



34.3.1 S'-matrix elements 

The energy correlation function of S-matrix elements is defined as follows 

Cf cd {uo) = (s$*(E x )Sg{E2j) = j™dte 2niult Cg cd (t) , (34.3.1) 

where Sg = S — S stands for the fluctuating part of the S-matrix. Being 
interested in local fluctuations (on the energy scale of the mean level spacing), 
it is natural to measure the energy separation in units of A. In the RMT limit 
N—>oo, function (|34.3.ip turns out to depend only on the frequency u = E2 ^ El . 
The Fourier transform Cg bcd (t) describes a gradual loss of correlations in time 
(Cg (t) = at t < identically, as required by causality). Physically, it 
is related to the total current through the surface surrounding the scattering 
centre, thus, a decay law of the open system |Lyu78 [DitOO]. 



The energy dependence becomes explicit if one considers the pole represen- 
tation of the S-matrix which follows from (|34.1.2jk S ab (E) = 5 ab - i 
Due to unitarity constraints imposed on S, the residues and poles develop non- 
trivial mutual correlations [Sok89] . For this reason the knowledge of only the 
JPDF of the resonances {£ n } (considered in the next section) is insufficient to 
calculate (|34.3.ip . The powerful supersymmetry method (see Chap. 7) turns 
out to be an appropriate tool to perform the statistical average in this case. 
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In their seminal paper |Ver85| . Verbaarschot, Weidenmiiller and Zirnbauer per- 
formed the exact calculation of (|34.3,ip at arbitrary transmission coefficients 
in the case of orthogonal symmetry. The corresponding expression for unitary 
symmetry was later given in [Fyo05|. These results are summarized below. 
The exact analytic expression for Cg bcd (uj) can be represented as follows: 



C abcd = s ab 5 cd T a T c ^/(l-T a )(l-T c )J ac + (5 ac 5 bd + S w S ad S bc )T a T b P ab . (34.3.2) 

(The same representation holds for Cg (t), so the argument can be omitted.) 
Here, the Sip term accounts trivially for the symmetry property S ab = S ba 
in the presence of time-reversal symmetry. Considering expression (|34.3.2|) in 
the energy domain, the functions J ac (^) and P a b(u) can be generally written 
as certain expectation values in the field theory (nonlinear zero-dimensional 
supersymmetric a- model). In the (3 = 1 case of orthogonal symmetry, one has 



/ / 2 \ / 2 \ \ §og 

Jac{u) = (IE 1+^7 + dfe) (E T+T^Jh + T^T^jFm) 

2 

p , (, ,\ — I ( V i 2 W (1- W ) \ TT,,\ g ° e 

Pab{u) ~ [ (^L (l+T oW )(l+T 6W ) + (l-T oWJ )(l-T 6/10 )J >M / M 

\2 



(34.3.3) 



where .F M = nj (i+rVo(i+r * s * ne so ~ cane d "channel factor", which 

accounts for system openness, and ((• • • ))^ oe is to be understood explicitly as 

1 f°°, f 1 d/xo(l-/io)A t o|/ii-/x 2 |e i7r ^ 1 +^+ 2 ^) 

8 7 ^Jo M Vo [(l + wVi(l + M2V 2 ] 1 / 2 (/io + /U 1 ) 2 (M + ^) 2l " j ' 

In the (3 = 2 case of unitary symmetry, the corresponding expressions read: 

gue 



Jac(u) = ( T7$— + 1-^— tttT + t-^— )T M 

\\ 1+Tafi! 1-TaW) ) \ 1+TcHi 1— TctM) J /a 



(34.3.4) 



j-m = n c fefe ((• • • ))r = ; «d M1 jo^ ^g^g^ c . .)• 

There are several important cases when the above general expressions can 
be simplified further. We mention first elastic scattering when only a single 
channel is open (see [Dit92] and Eq. (19) in |Fyo05|). The profile of the cor- 
relation function is strongly non-Lorenzian in this case. This is related to the 
power-law behavior of the form- factor at large timefl Cg bcd (t) ~ f-M/3/2-2 - n 
the general case of M open channels. Such a power-law time decay is a typi- 
cal characteristic of open chaotic systems [Lew91j lDit92j . Formally, it appears 



4 We conventionally measure the time in units of the Heisenberg time tn = 2n/ A {h = 1). 
The dimensionless time is t = fio + |(/ii + ^2) or t = /no + AH at /3 = 1 or /9 = 2, respectively. 
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due to the dependence Fm ~ FLU + § T c t)~~PI 2 of the channel factor at t^> 1. 
Physically, such a behavior can be related to resonance width fluctuations which 
become weaker as M grows, see Sec. 134.3.31 below. In the limiting case of the 
large number M > 1 of weakly open channels, T c <C 1, all the resonances ac- 
quire identical escape width Yl c =i T c , the so-called Weisskopf 's width (in units 
of so that Fm becomes equal to e~ t ^^ Tc . As a result, fluctuations in the 
S-matrix are essentially due to those in real energies of the closed counterpart 
of the scattering system [Ver86| . with the final expression in this limit being 

Cf cd (0j) = ( S ^ SM + S l^ ad5bC ) T aT b +S ab S cd TaTc r dt[1 _ b2{t)]e -(j:T c -2m^t 

2^T C - 2-kiuj J 

(34.3.5) 

where hi (t) is the canonical two- level RMT form-factor. 

As a straightforward application, we mention a prediction for the experi- 
mentally accessible elastic enhancement factor var(S aa )/var(S a ' ) ), a ^ b, which 
follows from the S-matrix correlation function at oj = 0. In particular, in the 
Ericsson's regime of many strongly overlapping resonances (^ T c S> 1) the first 
(dominating) term in (|34.3.5p yields the well-known Hauser-Feshbach relation. 
Very recently, the above results have been generalized to the whole crossover 
regime of gradually broken time-reversal symmetry (GOE-GUE crossover) and 
also tested in experiments with chaotic microwave billiards [Die09j . 

The above formulae can be also applied to half-scattering or half- collision 
processes (thus decaying quantum systems). Important examples of such events 
are represented by photo-dissociation or atomic autoionization where an absorp- 
tion of photon excites the quantum system into an energy region (with chaotic 
dynamics) which allows the subsequent decay. The optical theorem can be used 
to relate the autocorrelation function of the photo-dissociation cross-sections to 
that of S'-matrix elements (sec |Fyo98| for an alternative approach). Another 
particular interesting feature is the so-called Fano resonances resulting from an 
intricate interference of long-time chaotic decays with short-time direct escapes. 
This can be again described in terms of the S'-matrix correlations given above, 
see |Gor05| for explicit expressions, detailed discussion and further references. 



34.3.2 S'-matrix poles and residues 

Resonances. Within the resonance approximation considered the only singu- 
larities of the scattering matrix are its poles £ n = E n —^T n . Due to the unitarity 
constraint their complex conjugates £* serve as S-matrix's zeros. These two 
conditions yield the following general representation: 



(34.3.6) 
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which can be also verified using (|34.1.2|) by virtue of det S(E) = g^^r^fy- 
As a result, one can express the total scattering phaseshift (f)(E) = log det S(E) 
(or other quantities, e.g, time-delays discussed below) in terms of resonances. 
Thus statistics of the £ n 's underlay fluctuations in resonance scattering. 

To be able to analyse systematically statistical properties of N eigenvalues 
£ n of T~L e {i, it is natural to start with finding their JPDF Vm{£). Unfortunately, 
such a density is known in full generality only at M = 1. In this case the 
anti-Hermitian part of T~L e s has only one non-zero eigenvalue which we denote 
—in. Assuming k > to be non-random, one can find in the GOE case [Fyo99 

»*•-(«» * n hR ^rfH ■ (34 - 3 - 7) 

whereas in the GUE case the JPDF is given by 

^ 1 ({^})cxn^ Re ^|A({^})| 2 ^ r <j(« + ^Im^). (34.3.8) 

k=l ^ i=l ' 

Closely related formulae for randomly distributed k were first derived in [TJ1169J 
and, independently, in [Sok89| (see also Ref. [13] in |Fyo03| ). Actually, exploit- 
ing a version of the Itzykson-Zuber-Harish-Chandra integral allows one to write 
the analogue of (J34.3.8P at arbitrary M |Fyo99j |Fyo03| but the final expression 
is rather cumbersome. 

As typical for RMT problems, the main challenge is to extract the n-point 
correlation functions corresponding to Eq. (134. 3. 7j) or (I34.3.8p . The simplest, yet 
nontrivial statistical characteristics of the resonances is the one-point function, 
i.e. the mean density. Physically, it can be used to describe the distribution 
of the resonance widths T n in a window around some energy E. An important 
energy scale in such a window is the mean separation A between neighbouring 
resonances along the real axis. One can show [Som99] that Eq. (j34.3.7H implies 
the following probability density for the dimensionless widths y n = 7rr n /A: 

PmIM = / V - ^)e 2Xy F(X,y)d\ (34.3.9) 
with F(X,y) = (g-X) f°° % e -?ikM f? ^ ( P -r) eM-ry) where h cqu _ 

^ ' y U 9 (A-p) 2 V(p 2 -l)(p-g) J l (A-r)V(ra -l)(fl--r) 

pling constant g = 2/T — 1 > 1 is related to k by Eq. (|34.1.3|) . An important 
feature of this distribution is that it develops an algebraic tail oc y~ 2 in the 
case of perfect coupling g = 1. Such a behaviour was recently confirmed by 
experimental measurements of resonance widths in microwave cavities [Kuh08]. 

One can actually derive the distribution of the resonance widths for any 
finite number M of open channels. In the GOE case the corresponding expres- 
sion is rather cumbersome (see Eq. (3) in |Som99] ) but in the GUE case the 
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M-channel analogue of (|34.3.9|) is relatively simple, being given by |Fyo97a| 

for equivalent channels (all g c = g). This formula agrees with the results for 
chaotic wave scattering in graphs with broken time-reversal invariance [Kot00| . 
Note that in the limit of weak coupling g 3> 1 resonances are typically narrow 
and well-isolated, T n /A ~T<1. In this regime the above expressions reduce 
to the standard Porter-Thomas distributions derived long ago by the simple 
first-order perturbation theory (see discussion and Ref. [Por65] in Chap. 2). 

Actually, for the case of GUE symmetry and the fixed number M <C N — > oo 
of open channels, one can find in full generality not only the distribution of the 
resonance widths but all the n-point correlation functions of resonances S n in 
the complex plane. Assuming for simplicity a spectral window centered around 
Re£ = 0, all correlation functions acquire the familiar determinantal form 

lim^ -j^R n (z 1 =N£ 1 ,...,z n = N£ n ) = det [K(z h z*)]" fc=1 , (34.3.11) 

where the kernel is given by |Fyo99| 

M 



K(z 1 ,z* 2 ) = F 1 A { 2 (z 1 ) 



)Fll 2 (z* 2 ) d\e- iX ^-^H(g k + \), (34.3.12) 
J ~ l k=i 

with Fm(z) = Ylc=\ fl ~ (g -g C ) • ^ n ^ ne particular case of equivalent channels 
the diagonal K(z,z*) reproduces the mean probability density (|34.3.10p . 

Finally, we briefly describe the behaviour of S matrix poles in the semiclas- 
sical limit of many open channels M ~ — > oo. In such a case the poles form 
a dense cloud in the complex energy plane characterized by mean density p(z) 
inside the cloud [Haa92^ ILeh95a] . The cloud is separated from the real axis by 
a finite gap. The gap's width sets another important energy scale, determining 
the correlation length in scattering fluctuations [Leh95aJ ILeh9 5b| . Considering 
the M — > oo limit of Ki^zi^z^), one concludes that after appropriate rescaling 
statistics of resonances inside such a cloud are given by a Ginibre-like kernel 

\K{z l} z* 2 )\ =p(z)exp{-ivrp(z)|zi - z 2 \ 2 }, z = \{zi + z 2 ). (34.3.13) 

This result is expected to be universally valid for non-Hermitian random ma- 
trices in the regime of strong non-Hermiticity |Fyo03], see also Chap. 18. 
Nonorthogonal resonance states. As was mentioned in the introduction, 
the resonance states (quasimodes) are identified in our approach with the eigen- 
vectors of the non-Hermitian Hamiltonian % e ff- These eigenvectors form a bi- 
orthogonal rather than orthogonal set, the feature making them rather different 
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from their Hermitian counterparts. More precisely, to any complex resonance 
energy £ k correspond right \R k ) and left (L k \ eigenvectors 

n cS \R k ) = £ k \R k ) and (L k \n cS = £ k (L k \ (34.3.14) 

satisfying the conditions of bi-orthogonality, (L k \R m ) = 5 krn , and complete- 
ness, 1 = Ylk \ Rk){L k \. The non-orthogonality manifests itself via the matrix 
Omn = {L m \L n ) {R n \R m ) . In the context of reaction theory such a matrix is 
known as the Bell-Steinberger non-orthogonality matrix which, e.g., influences 
branching ratios of nuclear cross-sections [Sok89j and also features in the par- 
ticle escape from the scattering region |Sav97j . In the context of quantum 
optics O nn yields the enhancement (Petermann's or excess-noise factor) of the 
line width of a lasing mode in open resonators whereas O n ^ m describe cross- 
correlations between noise emitted in different eigenmodes (see Chap. 36). 

This attracted an essential interest in statistical properties of O nm which 
were first studied for Ginibre's complex Gaussian ensembles [Cha98| . For scat- 
tering ensembles, such a calculation at j3 = 1 can be successfully performed 
by re-summing the perturbation theory in the regime of isolated resonances at 
M = 1 [SchOO] or in the opposite case of many open channels, M S> 1, (strongly 
overlapping resonances) [MehOl]. The general non-perturbative result for the 
average density 0(£ ) = {J2 n ®n${£n — £)) of the diagonal parts is available at 
any M only for (3 = 2 symmetry [SchOO] ■ For the off-diagonal parts, the corre- 
lator £)(£,£') = {Yln^rn ®nm&{£n — £)$(£m — £')) is known non-perturbatively 
only for the single-channel case, see |Fyo02[ Fyo03| . 



34.3.3 Decay law and width fluctuations 

The decay law of open wave chaotic systems is intimately related to fluctua- 
tions of the resonance widths. This important fact can be best quantified by 
considering the simplest yet non-trivial decay function, namely, the leakage of 
the norm inside the open system |Sav97| (we use units of tjj = ^ = 1 below): 

P(t) ee l^Tr(e^«'e-^ff^ = -^(^CW^ - *" ') • (34.3.15) 

n,m 

In the case of the closed system P(t) = 1 identically at any time. Thus the entire 
time-dependence of the "norm-leakage" P{t) is due to the non-zero widths of 
the resonance states and their nonorthogonality. 

The exact analytic expression for Pit) obtained by supersymmetry calcula- 
tions turns out |Sav97"l ISav03| to be given by that for P a b(t), see Eqs. (|34,3.3|) 
and (|34.3.4|) . where one has to put T a = = appearing explicitly in the 
denominators. The typical behavior P(t) ~ Tl c (^ + ^^ct)~^^ 2 ■ In the so- 
called "diagonal approximation", which neglects nonorthogonality of the res- 
onance states and becomes asymptotically exact at large t, P(t) gets sim- 
ply related to the distribution of resonance widths by the Laplace transform: 
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fdiag(i) = ^(£„e" r ' li > = / °°dre- r V(r). With Eq. (134.3.101) in hand, one 
can then analyse in detail the time evolution in decaying chaotic systems and 
the characteristic time scales |Sav97j , see also |Cas97] for the relevant study. It 
is worth mentioning further applications of such a consideration to relaxation 
processes in open disordered conductors MirOO and nonlinear random media 
[Kot04| , to chaotic quantum decay in driven biased optical lattices [Glu02] and 
quasi-periodic structures ll\ol ()•"» . and also in the context of electromagnetic 
pulse propagation in disordered media [Cha03l ISki06| . 



34.4 Other characteristics and applications 
34.4.1 Time-delays 

Time delay of an almost monochromatic wave packet is conventionally described 
via the Wigner- Smith matrix Q(E) = -ihSUS/dE, see e.g. |Fyo97a| ldCn2] for 
introduction and historical references. This matrix appears in many applica- 
tions, e.g., describing charge response and charge fluctuations in chaotic cavities 
[Biit05]. In the considered resonance approximation of the energy-independent 
V the matrix Q(E) can be equivalently represented as follows (h = 1) |Sok97] : 

^ = v \e-h^e-h^ < 34 "> 

This gives the meaning to the matrix element Q a ^ as the overlap of the internal 
parts {E — / H c s) 1 V of the scattering wave function in the incident channels 
a and b. In particular, the diagonal entries Q cc can be further interpreted as 
mean time delays in the channel c. Their sum tw = jj^rQ can be related 
to the energy-derivative of the total scattering phase, tw = — jjjg log det S, 
and is known as the Wigner time delay. In view of (|34.3.6|) . this quantity is 
solely determined by the complex resonances: Tw(E) = jj £ n ^ E _ E T yi + ^ /a ■ 
Its average value is therefore given by the mean level density, (tw) = 2n/MA. 
Fluctuations of tw around this value can be quantified in terms of the au- 
tocorrelation function of the fluctuating parts, t w = tw — {tw)- The exact 
expression for f3 = 1 |Leh95b] and f3 = 2 |Fyo97a| cases are given by 

<&(E - §u)t*{E + _ f 2Re<(2 / , + m+/ , 2 ) 2 >J oe 

" { Re <( M0 + Mi) 2 )^ ' ( } 

(The same shorthand as in Eqs. (j34.3.3[) (j34.3.4[) has been used.) Actually, the 
corresponding result for the whole GOE-GUE crossover is also known |Fyo97b| . 

As it is clear from definition (|34.4.ip . the time-delay matrix at a given 
energy E is in essence a two-point object. Therefore, its distribution cannot be 
obtained by methods of Sec. 134.21 and requires the knowledge of the distribution 
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functional of S(E). The solution to this problem for the case of ideal coupling, 
S = 0, was given by Brouwer et al. |Bro97], who found the time-delay matrix Q 
to be statistically independent of S. Considering the eigenvalues q±, . . . , qu of 
Q (known as the proper time delays), the JPDF of their inverses j a = (^)qa 1 
was found to be given by the Laguerre ensemble as follows |Bro97| : 

M 

P(fh, . . . , 7m) cx J] \% -%fi[ t M,2 e-^l\ (34.4.3) 

a<b a=l 

The eigenvalue density can be computed by methods of orthogonal polynomials, 
explicit results being derived for M > 1 at any (3 and for f3 = 2 at any M. 
In the latter case, one can use (|34.4.3p for further extracting the distribution 
of Q aa and, finally, that of the Wigner time delay tw [SavOl] (this paper also 
provides those distributions at f3 = 1,4 and M = 1, 2). Representation (|34.4.ip 
together with the supersymmetry technique enables us to generalise the above 
results for the distribution of proper time-delays to the general case of non-ideal 
coupling [SomOlJ. In this case, the distribution of the so-called partial time- 
delays (related to Q cc ) is actually known in the whole GOE-GUE crossover 
|Fyo97b| . Finally, we note that in the special case M = 1 there exists a general 
relation between the statistics of the (unique in this case) time-delay and that 
of eigenfunction intensities in the closed counterpart of the system [Oss05j . The 
RMT predicts that such intensities are distributed for f3 = 1,2,4 according to 
the Porter-Thomas distributions, more complicated formulas being available for 
the /3 = 1 to /3 = 2 crossover |Som94|.[Fal94] . In this way one can easily recover, 
e.g., the above distribution of the time-delay in the crossover regime. The most 
powerful application of such a relation is however beyond the standard RMT for 
systems exhibiting Anderson localisation transition, see [Oss05t lMen05t [Mir06] . 

34.4.2 Quantum maps and sub-unitary random matrices 

The scattering approach can be easily adopted to treat open dynamical sys- 
tems with discrete time, i.e. open counterparts of the so-called area-preserving 
chaotic maps. The later are usually represented by unitary operators which 
act on Hilbert spaces of finite large dimension N, being often referred to as 
evolution, scattering or Floquet operators, depending on the given physical 
context. Their eigenvalues (eigenphases) consist of N points on the unit cir- 
cle and conform statistically quite accurately the results obtained for Dyson's 
circular ensembles, see e.g. [Ghi021 I()ss031 I.Iac031 ISchOQj for diverse models 
and physical applications. A general scattering approach framework for such 
systems was developed in [Eyo00| and we mention its gross features below. 

For a closed linear system characterized by a wavefunction ^ the "strobo- 
scope" dynamics amounts to a linear unitary map such that ^(ra + l) = u^(n). 
The unitary evolution operator u describes the inner state domain decoupled 
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both from input and output spaces. Then a coupling that makes the system 
open must convert the evolution operator u to a contractive operator A such 
that 1— A' A > 0. The equation + = A^(n) describes now an irreversible 
decay of any initial state \£(0) ^ when an input signal is absent. On the other 
hand, assuming a nonzero input and zero initial state ^(0) = 0, one can relate 
the (discrete) Fourier-transforms of the input and output signals at a frequency 
oj to each other by a M x M unitary scattering matrix S(u) as follows: 



S(lo) = Vl-rtf-f 1 " -ut, A = uy/l - fft , (34.4.4) 

e -iw _ a 

where f is a rectangular N x M matrix with M < N nonzero entries T%j = SijTj, 
< n < 1. This formula is a complete discrete-time analogue of Eq. (|34. 1.2j) . 
In particular, one can straightforwardly verify unitarity and show that 



det (At -e^) N ( Z *- P ^\ 

det S(cj) = e~ wN ) L = e -**N TT p. > (34.4.5) 

detfe---^ t^irt-thY 



where Zk stand for the complex eigenvalues of the matrix A (note \zk\ < 1)- 
This relation is an obvious analogue of Eq. (|34.3.6|) and gives another indication 
of Zk playing the role of resonances for the discrete time systems. 

Generic features of quantized maps with chaotic inner dynamics are emu- 
lated by choosing u from one of the Dyson circular ensembles. By averaging 
Eq. (|34,4.4p over u, one easily finds ftf = 1 — ((S 1 )! 2 . Therefore, M eigenvalues 
T a = t 2 < 1 of ftf play the familiar role of transmission coefficients. In the par- 
ticular case of all T a = 1 (ideal coupling), the non- vanishing eigenvalues of the 
matrix A coincide with those of a (N—M)x(N—M) subblock of it. Complex 
eigenvalues of such "truncations" of random unitary matrices were first studied 



analytically in ZycOO . The general ensemble of NxN random contractions 
A = uy/l — fft was studied in |Fyo03| starting from the following probability 
measure in the matrix space (cL4 = T\ dRe Aijdlm Aij): 

V{A)dA oc 5{A^A - G)dA , G = 1 - fft , (34.4.6) 

Equation (|34.4.6|) describes an ensemble of subunitary matrices A with given 
singular values. The question of characterizing the locus of complex eigenvalues 
for a general matrix with prescribed singular values is classical, and the recent 
paper |Wei08| provided a kind of statistical answer to that question. 

34.4.3 Microwave cavities at finite absorption 

The above consideration is restricted to the idealization neglecting absorption. 
The latter is almost invariably present, being usually seen as a dissipation of 
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incident power evolving exponentially in time. In the approximation of uniform 
absorption all the resonances acquire in addition to their escape widths T n one 
and the same absorption width r a b s = 7^ > 0. Operationally it is equivalent 
to a purely imaginary shift of the scattering energy E — ^ E + §r a b s = Ej. As a 
result, the S- matrix correlation function in presence of the absorption is related 
to that without by replacing u — > uo + i^j '2tt in (|34.3.1|) . In the time domain, 
the corresponding form-factor then simply acquires the additional decay factor 
e -7 * (in this way the absorption strength 7 can be extracted experimentally, 
e.g., for the microwave cavities [Sc h03] ) . 

Distribution functions undergo, however, more drastic changes at finite ab- 
sorption, as the scattering matrix S(E^) becomes subunitary. The mismatch 
between incoming and outgoing fluxes can be quantified by the reflection matrix 

R = S(E^S(E 7 ) = 1 - r abs Q(£ 7 ) . (34.4.7) 

This representation, with Q(E^) from (|34,4.ip . is valid at arbitrary T a b s |Sav03| . 
The so-called reflection eigenvalues of R play an important role for the descrip- 
tion of thermal emission from random media, as discovered by Beenakker, who 
also computed their distribution at perfect coupling (T = 1) (see Chap. 36.3). 
The above connection to the time-delay matrix at finite absorption allows one 
to generalize this result to the case of arbitrary coupling, T < 1 [Sav03| . 

In terms of quantum mechanics the quantity 1 — R aa has the meaning of 
probability of no return to the incident channel. For a single channel scattering 
without absorption it should be identically zero, but finite absorption induces 
a "unitary deficit" 1 — R aa > which can be used as a sensitive probe of the 
chaotic scattering. The quantity starts to play even more important role in a 
single-channel scattering from systems with intrinsic disorder. Indeed, for such 
systems the phenomenon of spontaneous breakdown of 5-matrix unitarity was 
suggested as one of the most general manifestations of the Anderson localisation 
transition of waves in random media, see |Fyo03b| for further discussion, and 
|Kot02| for somewhat related ideas. 

In the context of statistical electromagnetics, the matrix Z = iK(Ey) has 
the meaning of the normalized cavity impedance |Hem05]. In view of (134.1.11) . 
the joint distribution function of its real and imaginary parts fully determines 
the distribution of (non-unitary) S. It turns out that there exists a general rela- 
tion (fluctuation-dissipation like) between this function and the energy autocor- 
relation function of resolvents of K at zero absorption [SavQS] . This relation en- 
ables us further to compute the exact distributions of complex impedances and 
of reflection coefficients | S cc \ 2 at any M and for the whole GOE-GUE crossover 
Fyo05], and to explain the relevant experimental data [Kuh05b, Zhc06, Law08]. 
The corresponding study for the case of transmission is still an open problem. 

Quite often, however, an approximation of uniform absorption may break 
down, and one should take into account inhomogeneous (or localized-in-space) 
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losses which result in different broadening of different modes. The latter are 
easily incorporated in the model by treating them as if induced by additional 
scattering channels, see e.g. [ Roz03l ISav06bt IPol09| for relevant applications. 
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